Magnus Expansion and Three-Neutrino Oscillations in Matter * 
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We present a semi-analytical derivation of the survival probability of solar neutrinos in the three 
generation scheme, based on the Magnus approximation of the evolution operator of a three level 
system, and assuming a mass hierarchy among neutrino mass eigenstates. We have used an expo- 
nential profile for the solar electron density in our approximation. The different interesting density 
regions that appear throughout the propagation are analyzed. Finally, some comments on the 
allowed regions in the solar neutrino parameter space are addressed. 
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INTRODUCTION 



The need to introduce the three generations of neu- 
trinos into a computation of the transition probabilities 
of these particles while traversing a medium has been 
recognised long ago to try to accomodate the obser- 
vations of different experiments studying neutrino oscil- 
lations. Analytical treatments of three neutrino oscilla- 
tions in matter with varying density in the three gener- 
ation scheme have been studied in the past |2| aiming 
to deduce expressions for the oscillation probabilities of 
one type of neutrino into another. It has been shown, 
first in the case of two generations 0, Q| , and later in the 
case of three generations [j| that the differential equa- 
tion describing the evolution of the neutrino state in an 
exponentialy varying density profile could be solved ana- 
lytically in terms of confluent hypergeometric functions. 
Corrections to the mixing parameters in matter calcu- 
lated as series expansions have been performed by Fre- 
und 0, and a different perturbative analysis has been 
done by Narayan in pj . Global analyses of the recent ex- 
perimental dathave been extensively studied both, in the 
two and three generations cases [8] . In a previous paper 
by D'Olivo and Oteo @ an approximate expression to 
the evolution operator using the Magnus expansion was 
found, only for the non-adiabatic regime in the exponen- 



tially varying density profile. In this paper we present 
a complete semi-analytical computation of the evolution 
operator for neutrinos in the same density profile, using 
the Magnus expansion approximation which will work 
propperly in the case of adiabatic and non-adiabatic evo- 
lution of the neutrino state. This paper describes some 
important results of the M. C. Thesis work presented by 
Luis G. Cabral-Rosetti in Let H denote the Hamil- 
tonian of a quantum systen and U = U(t,to) the time 
evolution operator satisfying the Schrodinger equation 
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HU, U(to,h) 



(1) 



When H is independent of time, or more generally when 
[// dt'H(t'),H(t)} = 0, the solution of Eq.flJ is formally 

U = Exp[-i/h f* dt'H(t')}. Then it is natural to ask 
wether a solution of the form U = Expf2 would always 
be possible. A method for finding such a true exponential 
solution (without time ordering) is supplied by the Mag- 
nus Expansion (ME) • The magnus operator Q — \nU 
satisfies a differential equation which in turn is solved 
through a series expansion: Q = , fi n , where each 
term f2„ is of order hr n . The first two contributions are 
explicitly given by 



f dhH(h) , 
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° 2 = ~2ft2/ dtl / dt 2 [H{t 1 ),H(t 2 )] . (2) 
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Recursive methods to obtain the succesive terms have 
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approximate time-evolution operator obtained as U ps 
Uk = Exp(^Jj =1 f2„) will be unitary. Here we use the 
first Magnus approximant to obtain (approximate) ana- 
lytical solutions to the problem of 3 neutrinos oscillations 
in a medium with varying density as the Sun. 
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Figure 1: Behaviour of the matter mixing angles as functions 
of the effective potential V(t). 



II. SURVIVAL PROBABILITY FOR ELECTRON 
NEUTRINOS 

From the dependency of sin 20™ (t), and sin 20™ (i) 
on V(t) (see Ref. ^(J an d 113 for all details), we can 
distinguish five interesting regions where the oscillations 
behave differently. 

(1) Region of Extreemely Low Density 

If < V(t) < A 21 , both mixing angles, 9™(t), 
and 0™(t), are close to their values in vacuum (see 
Fig.©), this is 0™(i) -» 12 ; 0™(i) -> 13 , giving 



V a (t) 



and "P 7 it) 



(3) 



leading us to the well known result of vacuum os- 
cillations of three neutrinos 



(T{v e -> ^ e )) = c\ 2 c\ 3 + s\ 2 c\ 3 + s* 3 
(2) Low Density Resonance Region 



(4) 



For Vlt) ps A 21 , 0™(i) is still close to its 
vacuum value (6™ 3 (t) ~ 9 13 ), while 9™(t) is at the 
low density resonance, i.e. 9™ 2 (t) ~ making 



V a (t) 



and 



cos 2 IT(t) cos 2 6™ (t ) cos 2 0™(t o ) 

+sin 2 U(t) sin 2 0^(*o) cos 2 6™(t ) , (5) 



which leads to the expression 
1 1 



{V{v^vS) = 2 + 2 C 1 ^ 2 ^-) cos20 "i(to) cos20 12 , (7) 
where 

V, - V(t ) 



cos 29™ (t ) 



(8) 



and 



TV = sin 2 



(*E(fo) - «S(T))Exp (- «,)] ■ (9) 
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P 7 (t) -» sin 2 0™(t o ) , 



(6) 



Figure 2: The survival probability for an electron neutrino 
comming from the Sun vs. S m 21 /2E, with the parameters (a) 
sin 2 20i2 = 0.1, sin 2 20 13 = 0.01 y R = 10 3 ; (b) sin 2 20 12 = 
0.3, sin 2 20 13 = 0.1 y R= 10 3 . 

The quantity Vj represents the tramsition proba- 
bility between the states v 2m y v lm , and has been 
derived by D'Olivo Q in the case of oscillations 
between two neutrino species. The adiabatic result 
can be recovered by seting Vj — 0. For k, » 1, Vj 
is exponentially supressed as expected in the asimp- 
totic regime. On the other hand, for < 1 there 
are significant corrections to the adiabatic approxi- 
mation which reduce the effect of the resonant tran- 
sition. Eq.Q was first obtained by Parke for 
two species using the Landau-Zener approximation 
for the crossing probability: Vj — Exp ( — 
In the extreeme non adiabatic case, A 21 — > 0, mak- 
ing K t — > 0, and from Eq. JSJ it follows that 
^12(^0) — > § j making Vj — cos 2 12 . Introducing 
this value of Pj in Eq.JJJ) we recover the vacuum 
result for two neutrinos: 

{T{u.^u.)) = 1 - \ sin 2 20 l2 , (10) 

This should be contrasted with the result 

(P(u e -» v e )) = cos 2 12 , (11) 

predicted by the Landau-Zener formula (and, in 
general, by any result derived from the Dykhne's 



3 



formula), which deviates from the correct limit 
given by Eq. l|l(Jfl . when 8 12 is large. The correct 
value for the extreeme nonadiabatic case has been 
derived before, under the assumption that the tran- 
sition between the adiabatic eigenstates occurs in- 
stantaneously at the time t = t, [lfij . 



regions: A 31 > A 21 . 
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Figure 4: The la allowed region in the A21 vs sin 2 2#i2 plane 
corresponding to the range 0.25 < P e < 0.43 predicted by the 
combined results of SNO and Super Kamiokande. Notice the 
double MSW triangle structure. 



Figure 3: The survival probability for an electron neutrino 
comming from the Sun vs. Sm 21 /2E, with the parameters 

(a) sin 2 26 12 = 0.1, sin 2 20 13 = 0.001 and R = 100. (b) (4) High Density Resonance Region 
sin 2 26»i2 = 0.1, sin 2 2(9 13 = 0.3 and R = 10 3 . 



However, in this case the corresponding result 
for V c l approaches cos 2 8 12 as (^-§p) 2 instead of 
linearly, as in Eq. ©. Non adiabatic effects start 
to become important in this region when k, is 
comparable to 1, whenever the neutrinos cross the 

8 m 2 

low density resonance. If E < 2 v(tl) cos ^^12, 
Vj = 0, and the propagation of neutrinos will 
be adiabatic for n, < 1. For this reason, the 
asimptotic exponential expression (h — > 0) for 
Vc, must be modified by hand to consider this 
situation. An effective way of implementing such 
modification is to multiply by a step function 
&(V(t ) — A 21 cos 29 12 ), in such a way that the 
transition probability vanishes if a neutrino is 
produced after the resonance. It is worthnoting 
that such modification is not necessary with the 
Magnus result Eq. © given that, as a function of 

8 m 2 

™ 21 , the difference 0™ 2 (t ) — 0™(T) behaves as a 
continuous step. 



(3) Intermediate Density Region 

For A 21 <C V(t) <C A 31 , the mixing angles in 
matter are 0™ 2 (t) « |, and 6™(t) k 13 , making 



When V(t) w A 31 , the mixing angle 6*™(t) is 
approximately equal to while 0™ 3 (t) is at its high 
density resonance value of 0'" 3 (t) as j. In this case 



V a {t) -» cos 2 0™(t o ) cos 2 0™(t„) , (14) 



and 



Vj(t) -> cos 2 n(t) sin 2 6™ (t ) 
leading to 

1 1 



(15) 



(TV-*.)) = 2 + 2(1-2^ h ) cos 29? 3 (t ) cos 29 13 (16) 



where 

cos 20™ 3 (t o ) 



V h - V(t ) 



y/(v(t ) - v h y + bi 



(17) 



and 



(9? 3 (t )-9? 3 (T))Exp (-«,)] , (18) 



V a (t) ~ C 12 C 13 



(12) 



and 



V 7 (t) a s? 3 , (13) 

giving the result shown in Eq.(0}, provided there 
exists a clear separation between the resonance 



which again implies the two neutrino oscillations 
result. The crossing probability between the in- 
stantaneous eigenstates v 3m and v 2m given by Vj 1 
is perfectly analogue to that studied in the low den- 
sity resonance region. We can recover the adiabatic 
case if we set V c h — 0, and for n h S> 1, V c h is ex- 
ponentially supressed. The extreeme nonadiabatic 
case requires A 31 — 0, and from Eq. ijSJ) we have 
^(to) -> § leading to Vj 1 = cos 2 13 . All these ob- 
servations will lead us to an expression of the form 
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of Ea. l|10(l with 9 12 replaced by 8 13 . Non adiabatic 
effects become important when K h « 1 provided 
the neutrino crosses the High density resonance. If 
E < (cos 26 13 - Br 1 sin 2 6 l2 cos26> 13 ), with 

R = A 31 /A 2 i, no transitions between the instan- 
taneous eigenstates can occur (V c h = 0), and the 
propagation is adiabatic for K h < 1. Again, the 
difference 9™ 3 (t a ) — 6™(T) behaves as a continuous 
step giving Vj 1 the appropiate behaviour. 

















0.29 < P e 


<0.39 


, <E> 


= 10 MeV 


SNO+SK 




0.48 < P e 


<0.56 


,<E> 


= 1.1 MeV 


SAGE+GALLEX 




0.29 < P e 


<0.39 


, <E> 


= 5.5 MeV 


HOMESTAKE 















Figure 5: Overlap of the allowed regions for the experi- 
ments SAGE+GALLEX (0.46 < P < 0.66) and SNO+Super- 
Kamiokande (0.25 < P < 0.43). R = S = 250. 



Plots of the v e survival probability (Pr^—^)) as a 
function of are shown in Figs. 2-3, for different val- 
ues of sin 26 12 , sin 2 26> 13 , and R. We use the exponential 
profile H3 N e (r) = 245 Exp(- 10.54 r/R ) N Avo cm- 3 , 
where N Avo is the Avogadro's number, r is the radial dis- 
tance measured from the center of the Sun, and Rq is 
the solar radius (Rq — 6.96 x 10 5 Km). Except for those 
regions close to the center or the surface, this is a good 
approximation of the electron density in the Sun. We fur- 
ther assume that the v e are produced at r = 0.08635 Rq, 
where N e is the central electron density predicted by the 
Standard Solar Model (SSM). In Fig0|we show the al- 
lowed region in the A21 vs sin 2 2#i2 plane obtained by 
simply plotting the points in the plane wit survival prob- 
abilities lying within the la range of values extracted 
from recent joint analyses of the results of SNO and Su- 
per Kamiokande 18]. To produce this region we used 
the estimated value P e = (0.34 ± 0.05) for the survival 
probability, and an average neutrino energy of 5 MeV. 
We computed similar regions using the survival proba- 
bilities estimated by the experiments SAGE, GALEX 
and HOMESTAKE, taking the average energies for this 
experiments as in [l^ , and looked for an overlap of these 
regions in the parameter plane (see Fig. 5). No intention 
to give statistical significance to this region exists, but 
only it is shown that our result it is consistent with those 
achieved by rigurous analyses. 



(5) Extreemely High Density Region 



For V(t) 3> A 31 the (v e ) oscillations are strongly 
supressed due to the fact that 9^(t) 
0™ 3 (i) w §, making 



5, and 



V a (t) 



(19) 



and 



VJt) 



giving the result of Eq. J3J 



(V(l/ e ->■ V e )) = c\ 2 C 4 13 + s\ 2 c\ 3 



(20) 



(21) 
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